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ABSTRACT 

The NGC 5044 galaxy group is dominated by a luminous elliptical galaxy which is sur¬ 
rounded by ~ 160 dwarf satellites. The projected number density profile of this dwarf 
population deviates within ~ 1/3 of the virial radius from a projected NFW-profile, 
which is assumed to approximate the underlying total matter distribution. By means 
of a semi-analytic model we demonstrate that the interplay between gravitation, dy¬ 
namical friction and tidal mass loss and destruction can explain the observed number 
density profile. We use only two parameters in our models: the total to stellar mass 
fraction of the satellite halos and the disruption efficiency. The disruption efficiency is 
expressed by a minimum radius. If the tidal radius of a galaxy (halo) falls below this 
radius it is assumed to become unobservable. The preferred parameters are an initial 
total to stellar mass fraction of ~ 20 and a disruption radius of 4 kpc. In that model 
about 20 per cent of all the satellites are totally disrupted on their orbits within the 
group environment. Dynamical friction is less important in shaping the inner slope 
of the number density profile since the reduction in mass by tidal forces lowers the 
impact of the friction term. The main destruction mechanism is tide. In the preferred 
model the total B-band luminosity of all disrupted galaxies is about twice the observed 
luminosity of the central elliptical galaxy, indicating that a significant fraction of stars 
are scattered into the intragroup medium. Dwarf galaxy satellites closer to the centre 
of the NGC 5044 group may exhibit optical evidence of partial tidal disruption. If 
dynamical friction forces the satellite to merge with the central elliptical, the angular 
momentum of the satellite tends to be removed at the apo-centre passage. Afterwards 
the satellite drops radially towards the centre. 

Key words: methods: numerical - galaxies: clusters: general - galaxies: dwarf - dark 
matter. 


1 INTRODUCTION 

It is generally accepted that the galaxy distribution in rich 
clusters of galaxies closely follows the dark matter distri¬ 
bution. In other words the number density of galaxies ob¬ 
served in clusters is well descr ibed by a Navarro, Frenk 
and White profile (Na varro et alJll 997l). T his behaviour is 
established by observations (see e .g. Icarlberg et alJ Il997t 
iBiviano fc GirardU 120031 : |Lin_et_aP |2 004|) as well as theo- 
retical investigations (e.g ^Ga^^WrlTl^lKMl: iKr avtsov et alJ 
|2004 iDiaferio et alj|200il ; bpringel et alJl20oiF . In contra¬ 
diction to these outcomes the number density profile of 
the galaxies in the NGC 5044 group shows a flattening to¬ 
wards the centre, whi ch is not i n agreem ent with the stan¬ 
dard NFW profile llMathews et al.11200411 . The NGC 5044 
group o f galaxies resembles closely the properties of fossil 
groups llPonman et, alJll994l : lifones et all 1200(4 . 1200,4) . The 
central elliptical galaxy in NGC 5044, with luminosity Lb = 
4.5 x 10 10 I/Bf7) is surrounded by ~ 160 known dwarf sa tellites 
iFerguson fc Sandagdll99(l ICellone fc Buzzonill2005lf . At a 


distance of 33 Mpc (e.g. lTonrv et all20f)lil the absolute mag¬ 
nitudes of the non-central galaxies range from Mb = —19 to 
— 13, i.e. 2 x 10 7 < Z© < 9 x 10 9 Lbq- About 80 per cent of 
the group members appear to be spheroidal or dE galaxies. 

The basic differences between the dark matter substruc¬ 
ture (subhalo) distribution and the smoothly distributed 
component of cosmic matter assemblies is that substruc¬ 
tures lose mass due to tidal forces and are affected by dy¬ 
namical friction, whereas these two mechanisms have no 
impact on the smooth component. Both of these mecha¬ 
nisms can cause differences between the substructure distri¬ 
bution and the smoothly distributed matter component in 
groups and clusters of galaxies. Thus it is noteworthy that 
the galaxy number density profiles in rich clusters are in 
good agreement with the smooth density distribution. Dy¬ 
namical friction and tidal mass loss seem to balance each 
other. Both mechanisms acting at the same time result in 
a distribution of the substructures which resembles the un¬ 
derlying smooth distribution. However, the radial number 
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density profile of satellite galaxies in the NGC 5044 group 
does not follow the underlying matter distribution. We es¬ 
timate the expected radial distribution of galaxies in NGC 
5044 by following the orbits in the time dependent poten¬ 
tial well of the hosting halo. Since the masses of the dark 
halos that initially surround the stellar component of the 
satellites are not known well, we treat their total to stellar 
mass ratio (/* = Mt 0 t/M*) as a parameter in our computa¬ 
tions. A second parameter which goes into the integration of 
the orbits is the disruption radius Rdus■ Satellite halos that 
are truncated by tides down to this radius are assumed to 
reduce below the luminosity limit of the observations. We 
will show that the number density profile of the NGC 5044 
galaxy population can be explained by the action of dynam¬ 
ical friction and tidal forces if total disruption of galaxies is 
allowed. When tidal disr uption is consid e red it is no longer 
necessary to assume, as iMathews et al. 1200 411 have done, 
that the radial number density profile of satellite galaxies 
in NGC 5044 group requires an inherent deficiency of dwarf 
satellite galaxies at high redshifts. The tidally stripped stars 
may affect the intragroup mediu m, in particular t he abun¬ 
dance of heavy elements le.g. IZaritskv et alJl2004fl . 

This paper is organised as follows. § [ 5 ] focuses on the 
integration of the galaxy orbits within the time dependent 
potential well of the hosting halo. In § |3] the number den¬ 
sity profiles of the resulting satellite distribution at z = 0 
for different sets of parameters are discussed. § El investi¬ 
gates in detail the impact of dynamical friction versus tidal 
disruption. § El summarises the results obtained. 


2 DYNAMICAL MODEL FOR NGC 5044 


The model discus sed he re is a derivative of the 
model presented by IMathews et all J2004). Very recently 
IZentner et all ll2005ll demonstrated that this kind of semi- 
analytical model shows good agreement with state of the 
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lite galaxies is assumed to be proportional to that of the total 
gravitating matter, as if a fixed fraction of baryons formed 
into low-luminosity galaxies at a very early time. Thus the 
number of galaxies that have entered the virial radius of 
the host halo at a given time is proportional to the total 
amount of mass accreted at that time. Ho st and sat ellite 
halos are modelled according to findings of lWechsler et all 
(2002). A random B-band luminosity is assigned to every 
satellite galaxy in such a way, that their luminosity distri¬ 
bution matches the distribution observed in NGC 5044. Af¬ 
ter a galaxy has entered the virial radius of the NGC 5044 
group for the first time, its orbit is followed in the grow¬ 
ing potential well of the host. The integration scheme takes 
dynamical friction and tidal mass loss into account. Satel¬ 
lites are assumed to be lost if their tidal radius falls below 
a certain radius. For simplicity we assume that all stars are 
formed at redshift 2 * = 7 (we also employed z* = 3 and 11, 
which does not change the overall picture). Dark matter ha¬ 
los formed at that epoch have central densities comparable 
to that of the stellar densities observed. 


The underlying cosmology is the presently favoured con¬ 
cordance model with Q m = 0.3, Ha = 0.7, a Hubble con¬ 
stant of 70kms _1 Mpc _1 and a normalisation of the power 
spectrum of <rg = 0.9. 


2.1 Modelling the host halo 

Following the analysis of IWechsler et al] ( 2002 ’) the mass 
growth of the hosting halo can be described by 

M v {z) = M vfi e~ 2afZ (1) 


where a/ = 1/(1 + Zf ) defines the halo formation time. The 
expansion factor a/ at the time of formation for NGC 5044 
can be estimated by 


_ ciao 

C-vir — • 

a f 


( 2 ) 


If at present time ao — 1 a conc entration of c v j r = 11.1 
is assumed and ci is set to 5.125 llZentner et al.lf200 5) the 
expansion factor at formation turns out to be 0.46 which 
translates into a formation redshift of 2 / = 1.2. Knowing the 
evolution of mass (Eq.0 and concentration (Eq.[5J the virial 
radius at that redshift can be determined by the following 
formula. 


/ 3M(a) 1 \ 1/3 

V 47t A (a)p c (a) J 


( 3 ) 


where A(a) is a outcome of the spherical collapse approxi¬ 
mation fsee lEke et al.ll99~3 : lBrvan fc Normanll998ri . It gives 
the virialised overdensity with respect to the critical density 
pc(a) at that epoch. M(a ) is the time dependent mass spec¬ 
ified in Eq. 0 With the above recipe the dark matter profile 
and the gravitational potential of the host can be computed 
for any radius and redshift. The orbits of the satellite halos 
are computed within this potential well. 


2.2 Initial galaxy properties 

We assume that all satellite galaxies form at the same early 
redshift 2 * and that their space density is proportional to 
that of the dark matter. Thus the satellites enter the host 
halo at the same rate as the dark matter. With this simple 
model the host mass M Vi i when the ith satellite first entered 
the virial radius of the host can be found from 

Af v , i = (4:) 

where 0 < 7 Zm < 1 is a random number and M V} o = 3.9 x 

10 13 Mq is the current virial mass of the NGC 5044 group. 
The virial mass of the host at entry of the ith satellite can be 
translated into a redshift Zi via Eq. Qand a virial radius r v ,i 
by Eq. EH The velocity of the galaxy at the virial radius is 
the free-fall velocity from the turnaround radius rt = 2r v ,i. 

ttv,i = (GM v ,i/r v ,i) 1/2 (5) 

The stellar mass of the ith galaxy is 

m« }i = T B L B ,i ■ (6) 

Here the galactic luminosity is randomly selected from 

the power law part of the Schechter luminosity function 
dN/dLs oc Lg 1 ' 2 (which is satisfied in NGC 5044): 

L B ,i = [M 0 - 2 - n L {Li 0 - 2 - l 2 -°- 2 )] ” 5 , 


( 7 ) 
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with Li = 2 x 10 7 L_bq and L2 = 9 x 10 9 Lb© and 7 Zl is 
another random (0 < 7 Zm < 1)- For the stellar mass to light 
ratio Tb in Eq.|Hlwe adopt 

T fl = 0.71 x L 01 (8) 

as discussed hv lTrniillo et a.l] ( 2004 1. We assume that the 
galaxies form in isolation. Thus it is reasonable to embed 
them into a dark matter halo. Since the ratio between to¬ 
tal and stellar matter in isolated dwarf galaxies is not well 
determined, we introduce a parameter /* = Mtot/M* which 
allows us to vary the star formation efficiency. The total 
mass of the satellite including the dark matter halo is then 
given by 

mtot,i — Tb • (9) 

The density profile of all matter components - dark, stellar 
and gaseous - is assumed to follow a NFW -pro hle which is 
in agreement with recent studies bv IColi'n et ah ( 2004 *) who 
showed that the density profiles of dwarf (N-body) halos 
are well fitted by NFW-prohles. The concentration of each 
satellite halo accreting at redshift Zi is computed by means 
of Eq. |21 For this purpose ao has to be interpreted as the 
expansion factor at arrival (a 0 = 1/(1 + z,)) and a/ corre¬ 
sponds to the formation redshift of the halo a/ = 1/(1 + 2 *). 
As mentioned above, the choice of different formation red- 
shifts 2 * = 3,7,11, has marginal influence on the results. 
Subsequently, we assume z* = 7 for the global formation 
redshift of the satellite halos. This simplifications do not 
pay tribute to the dissipative baryonic processes, lik e rad i- 
ation or adiabatic contraction (see e.g. lGnedin et a 012004*1. 
Due to these processes the central parts of the satellite halos 
might be more tightly bound and the total tidal destruc¬ 
tion might become more difficult. This shortcoming can be 
compensated by the variation of the the second parame¬ 
ter introduced here, the disruption radius Rdis- If the tidal 
forces manage to strip the outer matter shells down to this 
radius, the satellite galaxy is assumed to vanish, i.e. fall be¬ 
low optical observability. We adjust this parameter to reach 
agreement with the global satellite distribution observed in 
NGC 5044, without tackling the complex processes close to 
the halo centre. 


2.3 Integration of the orbits 

We integrate the orbits of each satellite galaxy starting at 
the virial radius when the satellite first enters the virialised 
halo of the host. The absolute value of the initial velocity 
is computed under the assumption that the satellite has ex¬ 
perienced free-fall from twice the virial radius, Eq. [S] The 
absolute value of the velocity is distributed to a radial and 
a tangential component. The tangential comp onent is as- 
sign ed in accordance wit h the results bv lTormenl dl997l) (see 
also lZentner et, aI1l2005lh He gives a distribution of circular¬ 
ity e = J/J c of dark matter particles at the virial radius. 
Here J is the specific angular momentum of the dark mat¬ 
ter particle and J c is the maximum possible specific angular 
momentum of bound particles at that radius, t = 0 or 1 in¬ 
dicates radial or circu lar orbits, respective ly. Based on these 
results we model (see lMathews et alJl2004fl the probability 
density for circularity p' = dp/de with a third order polyno¬ 
mial, subject to the constraints that p^O) = p'(l) = 0 and 


normalised so the integral of p'(e) from 0 to 1 is unity. The 
resulting distribution is 

p(e) = oe 3 — (1.5a + 6)e 2 + (0.5a + 6)e (10) 

where 


a = 60(1 -2(e)). (11) 


The results of Tormen suggest that (e) ss 0.5, for which p'(t) 
becomes a symmetric quadratic. 

The orbits of the satellites are found by solving 


— = u and 
dt 


du G M(r) r / du\ 

dt r 2 r + V dt J df ’ [ ’ 


where M(r) is the NFW mass within the radius r of the 
host halo at that time. Orbits are computed from redshift 
Zi (redshift of entry) to z = 0. The deceleration by dynami- 
cal friction llChandrasekharlll 943t lOolni et a.lJIlOOfih can be 
described by 


= — u4-7tI n AG 2 mpu 3 [erf(X)- fj-Xe x 1 . 

7T 1 ' 2 

(13) 

Here p{r,t) is the local density of the host halo, u = |u| 
is the velocity of the satellite, X = u/(\/2a) and cr(r, t) is 
the mean velocity dispersion in NGC 5044. We assume that 
the total dispersion can be ap proximat ed by the cold dark 
matter dispersion pr esented bv lHoeft et, al.ld2004l) (see also 
MaHiew^tal 2004 1. We choose In A = 3 as suggested by 
Zhang et al (120021 which is the point mass approximation. 
This approximation is justified by the very efficient tidal 
mass stripping which effectively truncates the radius of the 
satellite halo. 

The tidal truncation radius is estimated by the Jac obi 
limit (see lBinnev fc Tremainelll987l : lHavashi et al ] na). 



rj 


TTltotji 

3A7(< D) 


1/3 

D 


(14) 


Here D is the distance of the satellite to the centre of the 
group, M(< D ) is the mass of the host within D and mtot,i 
is the current mass of the halo. If the resulting truncation 
radius r j is smaller than the current radius of the satellite 
halo, then the truncation radius is assumed to be the cur¬ 
rent halo radius and all the matter outside of rj is stripped 
instantaneously. In our model a minimum disruption radius 
Rdis is introduced. If the truncation radius falls below this 
radius the satellite’s luminosity is assumed to fall below ob¬ 
servability. We will show models for Rdis = 0, 2,4 and 6 kpc, 
respectively. For Rdis = 0 all the satellites survive. For com¬ 
parison the effective radius of the lowest luminosity dwarfs in 
NGC 5044 is approximately 4 kpc. We do not allow a change 
of the density or vel ocity profile of the satel lite halo fol¬ 
lowing the truncation. iKazantzidis et ah ( 12004 ) have shown 
that a cuspy density profile remains af ter mass is lost by tidal 
forces. However, tidal heating (see e.g. lGnedinl2003l) may re¬ 
duce the central density and therefore reduce the mass con¬ 
tained within the truncation rad i us. In fact by means of N- 
body simulations lHavashi et alJ (120041) and lKravtsov et alJ 
(200J) have shown that some mass is lost even within the 
tidal radi us of sub- h alos orbit ing in the potential well of a 
host halo. lHavashi et alJ ('’2 0041) compare different estimators 
for the truncation radius, concluding that the Jacoby limit 
most accurately estimates this additional tidal mass loss. 
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Figure 1. Comparison of the observed galaxy distribution in NGC 5044 (thick solid line) with the different model profiles. The various 
thin lines indicate alternate mass ratios between the total mass of the halo and the stellar component as indicated by the legend within 
each plot. The upper left-hand panel shows the resulting profiles if tidal mass loss is not taken into account. In the upper right panel 
tidal mass loss is tracked but the halos never are assumed to be disrupted. The lower panels depict the profiles if halo disruption is 
implemented. The disruption radii are assumed to be on left- and right-hand side 2 kpc and 6 kpc, respectively. 


Therefore we apply the Jacoby limit. Nevertheless, due to 
the dense baryonic cores, more mass may be retained than 
estimated from pure N-body investigations. We assume that 
the net trend of the different processes can be represented 
by an appropriate choice of the disruption radius Rdis that 
generates the final observed number distribution. Since the 
force exerted by dynamical friction is proportional to the 
mass squared, tidal mass loss strongly reduces the impact of 
dynamical friction. 


3 NUMBER DENSITY PROFILES 

In this section we present the resulting number density pro¬ 
files for different sets of parameters. As mentioned above 
we use only two parameters, firstly the initial total to stel¬ 
lar mass ratio /, to vary the contribution of the cold dark 
matter halo, secondly the disruption radius Rdis which mod¬ 
ulates the destruction rate due to tidal forces. Fig. 0 shows 
the results for four different tidal disruption/destruction ef¬ 
ficiencies. 

The upper left panel displays the radial distribu¬ 
tion resulting from dynamical friction without tidal mass 
loss . This scenario is rather unphysical since obse rvations 
(e.g lMaiewski et al ] l2004[jjrUidTa-P4 

S simulation s (e.g"lDiemanc HiMd . 20f^TKravtsov et alJ 

Nagai fc Kravtsa^l^OOfili have proven that tidal mass 
loss is significant. However it is interesting to note that the 
different masses of the satellite halos due to different /,,- 
values do not result in greatly different profiles. All profiles 
fall a little below the projected NFW-profile and are equally 
discrepant when compared to the observations. Dynamical 


friction shifts the whole distribution inwards without seri¬ 
ously affecting the overall shape. 

The upper panel on the right in Fig.Qshows the result¬ 
ing distribution if tidal mass loss is allowed but total tidal 
destruction is not allowed (Rdis = 0). This pre scriptio n is 
simila r to the treatment of the cluster galaxies in iGao et all 
1 20041 ). In agreement with their findings the galaxy num¬ 
ber density profile follows closely the NFW-profile (at least 
within the ~ | of the virial radius). However the result¬ 
ing profiles are not in agreement with observations of NGC 
5044. This difference emphasises that the shallow potential 
wells of the faint dwarf galaxies in NGC 5044 must be con¬ 
sidered. Observations addressing the dwarf populations in 
clusters of galaxies suggest that the dwarf ga laxies can be 
strong ly aff ected by tidal for ces. According to lDa.hlen et a.lJ 
i2004fl and iLin et akl ll2004l) the decreasing dwarf-to-giant 
ratio with increasing surface density in clusters of galaxies 
indicates that the high density environments are hostile to 
dwarf galaxies. 

The lower panels in Fig. 0 compare the the observed 
projected satellite distribution with model profiles for two 
different tidal destruction efficiencies, Rdis = 2 and Rdis = 6 
kpc. In the lower left-hand panel only the distribution with 
/* = 1 matches the observed distribution below 450 kpc. But 
the strong increase above 450 kpc and the fact that a mass 
to light ratio of 1 entirely neglects the contribution of cold 
dark matter gives evidence against this particular model. 
The model profiles in the lower right-hand panel agree with 
the observation for /* = 10, 20 and 40 below radii of 450 kpc, 
but obviously exceed the data and the NFW profile for large 
radii. (It should be mentioned that the observed distribu¬ 
tion is likely to be incomplete at large radii.) The census of 
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Figure 2. Comparison of the observed galaxy distribution in NGC 5044 (thick solid line) with the different model profiles for a disruption 
radius of 4 kpc. The various thin lines indicate alternate mass ratios between the total mass of the halo and the stellar component as 
indicated by the legend. The thick dashed line belonging to a stellar to total mass ratio of 20 marks the favoured model. 


surviving satellite galaxies in Tab. 0 indicates that a large 
fraction of satellites (~ 30 per cent) have been destroyed. 

Fig- m presents a model with disruption radius Rdis = 
4 kpc that matches the data for NGC 5044 very well. In¬ 
terestingly this Rdis turns out to be equal to the effective 
radius of the faintest galaxies observed in NGC 5044, if the 
minimum B-band luminosity Lb = 2 x 10 7 Lb© is trans- 
fo rmed into an effective r adius R e using the relation given 
in iDe Riicke et ah ( 12004 ). The total to stellar mass ratio 
which approximates the data best is /* = 20, or perhaps 
a bit lower since choosi ng f« = 10 results in very a similar 
behaviour. IPrada et alJ ; 20(135 have shown that the consid¬ 
eration of the dark matter halo of isolated galaxies leads 
to Mtot/Ls ~ 130 for galaxies within the luminosity range 
of 1 x 10 10 < I® < 4 x 1O 1O Z®0 corresponding roughly to 
/* = Mtot/M* = 16. Therefore the chosen parameter set 
seems to be in good agreement with recent observations. As 
indicated by the bold line in Tab.0 these parameters lead to 
a stellar loss rate of ~ 20 per cent. The total luminosity of all 
these lost galaxies is 9.4 x 10 10 I®g. This luminosity is twice 
as high as the luminosity of the central elliptical galaxy in 
NGC 5044, Lb = 4.5 x 10 10 Lb 0 , but only ~ 35 per cent of 
the total B-band luminosity in the group including all dwarf 
galaxies Lb, a ~ 13.7 x 10 10 Lbq (corrected for incomplete¬ 
ness, the corresponding B-band luminosity within r v i r in the 
best fitting model is Lb,a ~ 16.9 x 10 10 Lbq) and galax¬ 
ies disrupted in the model. Evidently, the m issing stars are 
scattered in the intragroup mediu m (see e.g. lZaritskv et all 
120041 : ISommer-Larsen et al ]Eoq3). 

We have shown the results for a uniform formation red- 
shift of 2 * = 7. Different formation redshifts result in dif¬ 
ferent concentrations of the satellite halos and thus might 
affect the tidal stripping efficiency. Therefore we have re¬ 
done the modelling with formation redshifts of 2 * = 3 and 
11. The differences are marginal. We conclude this section by 


Rdis 

Mtot/M * 

-N cen / -N v i r 

Nt ■ 

1 y jnc 

Ntid 

Lbq 
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1 

0.26 

0 

173 

6.4 x 10 10 

2 

10 

0.10 

0 

63 

3.8 x 10 10 

2 

20 

0.08 

2 

47 

4.9 x 10 10 

2 

40 

0.06 

8 

31 

6.5 x 10 10 

4 

1 

0.68 

0 

486 

2.4 x 10 11 

4 

10 

0.25 

0 

170 

8.4 x 10 10 

4 

20 

0.21 

3 

137 

9.4 X 10 10 

4 

40 

0.18 

7 

110 

1.1 x 10 11 

6 

1 

1.34 

0 

1075 

6.3 x 10 11 

6 

10 

0.42 

0 

297 

1.5 x 10 11 

6 

20 

0.32 

3 

219 

1.3 x 10 11 

6 

40 

0.28 

7 

179 

1.4 x 10 11 


Table 1. Survival rates. The first two columns display the dis¬ 
ruption radius and total to stellar mass, respectively. The third 
column gives the fraction of galaxies that ended in the centre to 
the current number of galaxies within the virial radius. The fol¬ 
lowing two columns list the number of galaxies that disappeared 
due to dynamical friction and tidal disruption. The last column 
shows the summed luminosity of all galaxies that have been de¬ 
stroyed. The bold face line indicates the favoured model. 

repeating that the model parameters /* = 20, Rdis = 4 kpc 
and 2 * = 7 match the data very well. 


4 DYNAMICAL FRICTION AND TIDAL 
MASS LOSS 

In this section we show the orbital behaviour of three par¬ 
ticular galaxies. These three cases are chosen to make some 
of the features discussed above more explicit. In Figs. emu 
and 0 we compare the behaviour of satellite halos having 
the same initial conditions with (solid and dashed lines) and 

























6 


A. Faltenbacher and W.G. Mathews 


1.00 



1 10 100 
Distance [kpc] 


Figure 3. Trajectory of a model galaxy. The ratio of forces 
caused by dynamical friction and by gravitation is ploted versus 
the physical distance of the galaxy from the centre of the group. 
The large, right dot marks the first penetration of the galaxy 
through the virial radius. The dotted line represents the path of 
the galaxy if tidal mass loss is suppressed but retaining dynamical 
friction. The solid and dashed (almost totally hidden) lines show 
the path if tidal mass loss is also taken into account. The left 
dot marks the location where the galaxy is assumed to be tidally 
lost if a disruption radius of 4 kpc is applied. The dashed line 
follows the path of the galaxy till z = 0 if the galaxy suffers 
no tidal destruction (Rdis = 0). The solid horizontal line at the 
top separates the friction dominated area from the area where 
gravitation prevails. 

without (dotted lines) tidal mass loss. The dashed lines indi¬ 
cate the orbits without tidal destruction, whereas the solid 
line ends (additionally marked by a second large dot) when 
the truncation radius falls below 4 kpc. In these figures the 
ratio of the forces due to dynamical friction and gravitation 
( F/G ) is plotted versus the distance of the satellite to the 
centre of the group. The horizontal line for F/G = 1 indi¬ 
cates the threshold where dynamical friction equals the force 
exerted by gravitation. Once an orbit transgresses above this 
line, it is dominated by dynamical friction and the galaxy 
drops nearly radially toward the centre. 

The large, right dot in Fig. |3] marks the first passage 
through the virial radius and therefore the beginning of in¬ 
tegration. After orbiting once around the group centre the 
satellite moves farther outwards bacause its initial velocity 
was the free-fall velocity from twice the virial radius. The 
eccentricity of the orbit without tidal stripping (dotted line) 
changes from initially t ~ 0.6 to a final value of « 0.7. Inter¬ 
estingly the contribution of friction to the total force shows 
a minimum for intermediate radii and rises again if the satel¬ 
lite moves further out. Dynamical friction is not negligible 
during the apo-centre passage. For small velocities F/G is 
roughly proportional to r/u, where r is the distance to the 
group centre and u is absolute values of the velocity. Fric¬ 
tional decay of the tangential velocity at apo-centre leads 
to more radial orbits, if tidal mass loss is not taken into ac¬ 
count. Furthermore the computed orbit indicates that satel¬ 
lites can first transgress the F/G = 1 line during their apo- 
centre passage. In that case all the angular momentum of 
the satellite is lost and it approaches the centre on a pure 
radial orbit. Such an orbit might be appropriate for a mas¬ 



Distance [kpc] 

Figure 4. Trajectory of a model galaxy, the description is the 
same as in Fig. |2] The absence of the second black dot indicates 
that this halo is not destroyed by tidal forces when a disruption 
radius of 4 kpc is applied. 


sive black hole that experiences dynamical friction without 
tidal losses. The eccentricity of the orbit with tidal stripping 
(solid/dashed line) stays « 0.6 all the time. This is due to 
that fact that tidal stripping reduces the action of dynam¬ 
ical friction. Most of the mass is lost during the first infall, 
as can be inferred by the shallow slope till the first peri- 
centre passage. The large, left dot indicates where the tidal 
destruction is assumed by our approach. Here Rdi s = 4 kpc 
was chosen. Fig. [3 shows that after ~ 10 orbits the satellite 
is destroyed at a distance ~ 25 kpc from the centre, which 
is quite close to the central elliptical galaxy. 

The orbit shown in Fig. 2] has a similar behaviour, but 
differs from the previous orbit in that the satellite is not de¬ 
stroyed by tidal forces (no second large dot). Without tidal 
stripping and mass loss (dotted line) the satellite halo even¬ 
tually falls toward the group centre as indicated by the final 
passage above the F/G = 1 line. This catastrophic loss of 
angular momentum by dynamical friction occurs during an 
apo-centre passage. When all its orbital angular momentum 
is lost, the satellite approaches the centre on a radial orbit 
starting at ~ 60 kpc from the centre. The eccentricity of the 
orbit without tidal stripping changes from ~ 0.8 to « 0.7, 
whereas the orbit with tidal mass loss retains its initially 
eccentricity of « 0.8. 

Fig. 0 is an example for a very radial orbit in which the 
eccentricity remains at « 0.9 for all orbits, either with or 
without tidal mass loss. The large initial radius ~ 450 kpc in¬ 
dicates that the satellite approaches rather late, that means 
the potential well of the hosting group is already very deep. 
The large, left dot marks the destruction during the first 
peri-centre passage, if tidal stripping and destruction are 
allowed (Rdis = 4kpc). For this particular case, which is 
rather extreme, the fast destruction seems reliable, since the 
satellite is falling from ~ 900 kpc (twice the virial radius at 
epoch of entry) down to ~ 20 kpc which is very close to the 
central elliptical galaxy. 

The detailed analysis of the individual orbits clearly 
demonstrates that the influence of dynamical friction is 
largely reduced by tidal mass loss. The estimation of the 
tidal radius by the Jacobi limit and the instantaneous strip- 
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Figure 5. Trajectory of a model galaxy, the description is the 
same as in Fig. [21 

ping up to this radius lead to a rapid decrease of the halo 
mass during the first infall. The instantaneous stripping de¬ 
scription might reduce the mass a bit to fast, but after very 
few orbits the Jacobi radius should be a good approxima¬ 
tion. On subsequent orbits dynamical friction tends to push 
the halo closer to the centre which is accompanied by fur¬ 
ther tidal mass loss, since the tidal radius of the halo is 
proportional to its distance from the centre. The force ex¬ 
erted by dynamical friction is proportional to the halo mass 
squared, thus the loss of mass reduces the impact of friction 
further. In conclusion the dominant process for shaping the 
inner number density profile of groups like NGC 5044 is tidal 
disruption and/or destruction. 


5 CONCLUSION 

Using a simple model we find that the interplay between dy¬ 
namical friction and tidal disruption can explain the shallow 
slope of the projected number density profiles of satellite 
galaxies of the fossil group NGC 5044. The strong deviation 
of observed satellite galaxies from NFW profiles is a result 
of the tidal disruption of the stellar cores of galaxies orbit¬ 
ing closer to the group centre. This feature may be more 
naturally explained in this manner than by assum ing an 
abs ence of dwarf galaxies at high redshifts ((Mathews et alJ 
l2004lf . We use two parameters to model the observed dwarf 
galaxy distribution: the total to stellar mass ratio /* and 
the disruption radius Rdis■ With parameters /* « 20 and 
Rdis ~ 4 kpc the model can match the observed projected 
radial number distribution in NGC 5044 very well. J* m 20 i s 
in good agreement with the findings of lPrada et alJ (2003)• 
They show that the Mtot/Ls of isolated galaxies with 13- 
band luminosities of 1 x 10 10 < Lb < 4 x 10 10 Z®@ is in 
the range between 100 and 150. For /* = 20 we obtain 
Mt ot/Ls ~ 160. If we us e the R e — Lb relation presented 
by EE Riicke et alJ i 2 1111 l.i to compute the effective radius 
for our lowest luminosity model galaxies (Lb ~ 2x 10 7 Lbq) 
we obtain R e « 4 kpc. Therefore, the global galaxy distribu¬ 
tion in NGC 5044 seems to favour a destruction radius that 
is close to the minimal R e of the group satellites. 

We find that the total luminosity of the tidally dis¬ 


rupted galaxies is about twice the luminosity of the cen¬ 
tral elliptical or ~ 35 per cent of the total optical light of 
the group, allowing for incompleteness of identified satel¬ 
lite galaxies. This indicates that a large fraction of satel¬ 
lite stars are scattered into the intragroup medium by 
tidal fo rces. Using t he near-infrared A'-band properties of 
clusters iLin fc Mohr : (200 4) estimate that a large fraction 
(50 per cent) of the total stellar luminosity is contributed by 
i ntracluster light. This fraction increases with cluster mass. 
iNeill et alJ (2005) searched the Fornax cluster for novae in 
between the galaxies and find that ~ 16 — 41 per cent of the 
total cluster light comes from the intracluster stars. Dwarf 
galaxies orbiting close to the centre of groups like NGC 5044 
may show morphological evidence of partial tidal destruc¬ 
tion. 
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